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Abstract 

We give an elementary derivation of the chiral magnetic effect based on a strong 
magnetic field lowest-Landau- level projection in conjunction with the well-known axial 
anomalies in two- and four- dimensional space-time. The argument is general, based on 
a Schur decomposition of the Dirac operator. In the dimensionally reduced theory, the 
chiral magnetic effect is directly related to the relativistic form of the Peierls instability, 
leading to a spiral form of the condensate, the chiral magnetic spiral. We then discuss 
the competition between spin projection, due to a strong magnetic field, and chirality 
projection, due to an instanton, for light fermions in QCD and QED. The resulting 
asymmetric distortion of the zero modes and near-zero modes is another aspect of the 
chiral magnetic effect. 
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1 Dirac Operators, Dimensional Reduction and Axial 
Anomalies 



In this first section we present an elementary relation between the chiral magnetic effect 
[H El El m El E] and the axial anomalies in four-dimensional and two-dimensional space-time 
[Zl [8]. This follows from the basic structure of the Dirac operator, together with the lowest- 
Landau level (LLL) projection produced by a strong magnetic field, a basic feature of the 
phenomenon of magnetic catalysis [9j. We first review the sub-block structure of the Dirac 
operator and the associated Schur decomposition [TU] of the propagator. 



1.1 Lowest Landau Level Projection 

We adopt the following conventions for Dirac matrices in four dimensional Minkowski space: 

where a-', j = 1,2,3, are the 2x2 Pauli matrices, and 1 is the 2x2 identity matrix. It is 
convenient to write these 4x4 Dirac matrices in 2 x 2 block form: 

^-(;:) (^) 

where the 2x2 matrices and are: 

a^" = (1, -a^) , = (1, a^) (3) 

The Dirac matrices satisfy the anti-commutation relations, {7^,7^^} = 2//^^, with Minkowski 
metric rj^'^ = diag(l, —1, —1, —1). 

The Dirac operator is defined tohe = 7^^ (9^ — ieA^j), where later we will take the gauge 
field to have both abelian and non-abelian parts, but for now we take it to be abelian. 
Thus, we have a natural decomposition of the 4x4 Dirac operator into 2x2 sub-blocks: 



V = -iO + m 



m D 
-D m 



(4) 



Since det {—il^ + m) = det {il^ + m), we consider 



det (z + m) det (-z ^ + m) = + ^ ^ (5) 

The operators D D and D D differ from one another, and from the Klein-Gordon operator 
Dfj^D^, in their spin-projection terms: 

DD = D^D^^+^-d^^F^, (6) 

DD = D^D^ + ^a^^'F^, (7) 



1 



where the 2x2 spin matrices, 

a^*^ = 1 {a^a" - a^a^") , c^''" = ^ (a^a" - a'a^") (8) 

are the sub-block components of the usual 4x4 spin matrices 

1 /cr^'' \ 



Now, suppose we have a background electromagnetic field consisting of a magnetic field 
B = (0,0, -B), and an electric field E = (0,0, ii^), both of which are pointing in the 
direction. We do not need to assume that these fields are uniform, but we choose B = 
B{x'^,x'^) and E = E{x°,x^), choosing gauge field = {A^ (x^) , {x'^) , {x'^) , (x^)) , 
which satisfies duA^^ = 0. Then, with only a third component of both the electric and 
magnetic field, the 2x2 operators in ^ and ([T]) reduce to 

DD = D^D^" - e{B + iE)a^ (10) 
DD = D^,D^ - e{B -iE)a^ (11) 



because = = —cr^, while a^^ = —a^^ = ia^. The important point in (10) and (11) 
is that both 2x2 operators D D and D D are diagonal in terms of their Dirac matrix 
structure. 

To understand the effect of a strong magnetic field, and in particular its projection to the 
lowest Landau level (LLL), consider the factorization of the Klein-Gordon operator, for this 
case parallel E and B in the x^ direction: 

D^D^ = {Do TD;){Do±D^)±ieE- {D, t ^^^2) (D, ± 2^2) ±eB (12) 

If the magnetic field is constant, and for example in the symmetric gauge Ai = —^x"^, and 
A2 = §x^, then we can adopt complex coordinates, z = x^ + ix"^, z = x^ — ix"^, to write 



D,±,D, = rsrj„: (13) 




If eB > we choose the upper sign to obtain normalizable solutions, with the Gaussian 
factor exp(— ^|zp), in which case we see that the eB term in ( [l2| cancels the spin term in 
(10) and (11), when = +1. Thus the magnetic component of the Dirac operators leads 
to a zero mode of the 2x2 operators DD and DD, when the spin is aligned along the 
direction of the magnetic field. In fact, this also applies to the situation of an inhomogeneous 
i?(a;^,a;^) field, and the degeneracy of this lowest-Landau-level (LLL) is given by the integer 
part of the net magnetic fiux; this is the Aharonov-Casher theorem pT] [12]. the projection 
onto the LLL. 

1.2 Schur decomposition of Dirac propagator 

The Schur decomposition gives an elementary algebraic decomposition of the inverse of a 
matrix in terms of its sub-block structure [TU]. This leads immediately to an associated 
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sub-block decomposition of the Dirac propagator, and as we show, it also provides a simple 
description of the lowest Landau level projection in a strong magnetic field, which is the key 
to magnetic catalysis [9]. 

Consider a matrix M written as 

where a and d are square, but not necessarily of the same size [and so h and c need not be 
square matrices], then we can write the inverse of M in block form in two different ways: 

= i-d'^cs-^ d-^ c7-Hd^^ + d-^) ' '^^-bd-'c (15) 



t = d-ca-^b (16) 



where s and t are the two different Schur complements of M. Note that the first expression 
only requires d and s to be invertible, while the second expression only requires a and t to 
be invertible. Applying this to the Dirac operator in Q, we find 

s= — (m^ + DD) , t = — (m^ + DD) (17) 
m \ / m \ / 

which gives then two different decompositions of the Dirac propagator: 



l-D ^\^ d) D 



1 jj 



(19) 

The Euclidean analogue of this decomposition corresponds precisely to the chiral decompo- 
sitions used in [131 E] • 

(20) 

where Q^^"* = l/(m^ + D D), and Q^^^ = l/(m^ + D D). These chiral decompositions of the 
fermion propagator are particularly useful since they show clearly the projection onto chiral 
zero modes. 

In order to use such Schur decompositions to characterize also the lowest Landau level 
projection, in addition to the chiral decomposition, it is convenient to write the propagators 



in (18) and (19) in factored form: 



= - ( 1 n J ( I) ( ; 1 ) (21) 



1 







m 


\ m 


1 


1 


il 


m 


m 




1 



0\ /I -^D 



1 \ / 1 



-^Stt^i y-^D 1 



(22) 



We stress that at this point wc have only used elementary algebra to express a general Dirac 
propagator in terms of its 2 x 2 sub-block structure. 

From these Schur decompositions, we can reduce the expectation values of the charge 
and axial currents to much simpler forms. By straightforward manipulations we find for the 
charge current 

<f> - ^etr4x4(7''2^"') (23) 
= i etr2x2 (( Da''- a'' d)—^—^] (24) 



= 2e?7'^"tr2x2 (^A 



1 



m? + DiD, 



(25) 



and for the axial current 

<ir> = ietr(7^75P-^) (26) 
= I e tr2x2 i{Da'' + a'^D) ^- S) (27) 

V + D D / 

Here we have used the facts that {a" a'* + a'' a") = 2?]'"' l2x2, while {a" a'* - a'^ a'^) = 2ia'"'. 
In particular, note that since a^^ — —ia^ — — cr^°, 

<f> = 2e77^'-tr2x2fl^o (29) 

<f> = -2e77^-tr2x2 fi^3 , (30) 

V -\- D D J 

<jI> = 2er;^-tr2x2fa^i^3 , ,\ ^ ) (31) 

<jI> = -2e7^^^tr2x2f<7'l^o— ^--') (32) 

V + D D / 

Again, thus far we have only used elementary algebra to reduce the expectation values of the 
4x4 matrices and propagators to expressions involving just 2x2 matrices and propagators. 



1.3 Currents and anomalies in the lowest Landau level projection 

Now suppose the background field consists of a very strong magnetic field in the direc- 
tion. Then we project onto the LLL, which means that the states contributing to the 2x2 
expectation values have a"^ = +1. Therefore, we see immediately that in this LLL projection 
limit, where we project onto motion along the magnetic field direction, we can write for the 
remaining currents: 

<#>=e''^<Jiv> , M, TV = 0,3 (33) 
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where the epsilon symbol is = +1 = — e^*^. This is exactly the relation between the charge 
and axial current in two dimensional space-time. 

Furthermore, suppose the four dimensional background field consists of a strong magnetic 
field and also an electric field, both directed along the axis. Then, following the analysis 
of the first section, we choose a gauge field of the form A^^ = {A^{x^), y4^(x^), A'^{x^), A'^{x^)), 
satisfying df^A'^ = 0. Then a simple two dimensional computation yields current expectation 
values for Dirac indices and 3 (we adopt the convention of using capital Roman indices 
M,N to denote the components of the dimensionally reduced {x^,x^) plane): 



where |^ is the Landau degeneracy factor, in its local Aharonov-Casher form. Note that 
these expressions are consistent with charge current conservation and the two dimensional 
axial anomaly: 

9m<J^>=0 , dM<J^>='^— (35) 

In this LLL projection limit, we can alternatively express this result in four dimensional 
language as 

2 2 

d, < f >LLL= , d^< >LLL= ■ E = ^F'^'F^. (36) 

which expresses charge current conservation and the four dimensional axial anomaly. This 
makes it clear that the relevant anomaly is the "covariant" anomaly, rather than the "con- 
sistent" anomaly. For abelian theories these differ by a factor of 1/d in {2d — 2) space-time 
dimensions, while for non-abelian theories the covariant and consistent anomalies have dif- 
ferent field structure [151 [IB] • 

To make the connection with the chiral magnetic effect, we note that it is natural to iden- 
tify A'^{x^) with a spatially dependent chemical potential, and A^{x^) with a time dependent 
chiral chemical potential: 

A° ^ /i , o /is (37) 

For A^ and /i this is obvious, because the coupling is given by fiip'j^ip- For A^ and /i5, this 
follows because the coupling is ^^ip'j^'j'^ip: and in the LLL projection 757° -H- 7^, since 7^ 
has off-diagonal sub- blocks =FC"^, and cr^ — >■ +1 in the LLL limit. 



Therefore, we can understand the two dimensional currents in (34) as 



n o U5 eB , . 

<^I>='tiI ' <A>=^^ (39) 

These relations express the chiral magnetic effect, which we see is a direct consequence of the 
axial anomalies in two and four dimensional space-time, after the LLL projection caused by 
a strong magnetic field. The coefficients are fixed completely by the anomaly equations. For 
a complementary discussion of the relation between chiral asymmetry and the axial anomaly, 
see di]. 
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Figure 1: Parallel electric and magnetic fields produce electron-positron pairs from vacuum, 
and because of the LLL projection caused by the strong magnetic field, correlated with spin 
and charge, this results in a net flow of chirality along the direction of the fields, in accordance 
with the Schwinger pair production rate and the chiral magnetic effect. 

1.4 Chiral Magnetic Effect and the Schwinger Effect 

The chiral magnetic effect can also be understood naturally in terms of the Schwinger effect 
[T8| [19], particle production from vacuum, which occurs when there is a non-zero electric 
field background, as illustrated in Figure [TJ With approximately constant parallel electric 
and magnetic fields directed along the axis, the Schwinger pair production rate, per unit 
volume, is given by 

r = e^||cothf|.)e-">VI-l (40) 



47r2 \E 

When i? — 7- we recover the usual result for a pure electric field, and for nonzero magnetic 
field we find an enhancement of the rate which is linear in B in the strong magnetic field 
limit. Positrons are accelerated along the direction of the electric field, and electrons in the 
opposite direction. In the massless fermion limit this corresponds to a net production and 
flow of chirahty, because the lowest-Landau-level projection projects spin according to the 
charge and the direction of the magnetic field. Thus we find the rate of change of chirality 

f = 2r = ^ (41) 

dt 27r2 ^ ' 



in agreement with the axial anomaly (36) and the chiral magnetic effect (39). The electric 
field produces the acceleration while the strong magnetic field provides the LLL projection 
that correlates spin with the direction of flow of charge, and hence also of chirality. Physically, 
a spatially dependent A'^{x^) produces charge separation, as for a local chemical potential, 
while a time dependent y4^(x°) drives a current along the direction of the electric field [20] . 
which in the LLL projection corresponds to a flow of chirality. 

1.5 Maxwell-Chern-Simons theory and the Schwinger model 

The interpretation of the chiral magnetic effect in terms of the effect of an electric field 
directed along the same direction as the strong magnetic field is also very natural in terms 
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of an effective Maxwell-Chern-Simons theory resulting from an adiabatic space- or time- 
dependent theta parameter p]. Express the theta term in the Lagrangian as (up to a total 
derivative) 

-^,dF^'^F^^ = P,rcs (42) 

where 

P, = d,e , J^^ = |1 e'^-P-A.F,. (43) 

The pseudo vector encodes the anomalous terms from the chiral magnetic effect, modifying 
the usual inhomogeneous Maxwell equations to read 

d^F^^- = r -^^P^F^^- (44) 




Figure 2: The effect of a spatially inhomogeneous theta parameter, with the interpretation 
of its gradient being a spatially inhomogeneous chemical potential /i ~ The resulting 
spatially inhomogeneous electric field produces a build up of negative charge at the left-hand 
inhomogeneity and positive charge at the right-hand inhomogeneity, producing the electric 
charge separation of the chiral magnetic effect. 
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A^{x^) = doe 



yU5(x°) 



Figure 3: The effect of a time dependent theta parameter, with the interpretation of its 
gradient being a time-dependent chiral chemical potential /is ~ in the LLL projection. 
Via the Schwinger effect, the resulting time dependent electric field produces a flow of electric 
current linear in A'^{t), and in the LLL projection this corresponds to a net flow of chirality. 



From (39), in the strong magnetic field limit we project to the 2d {x^,x^) plane with the 
identifications 



A, 



P. = An 



pM ^ j\ 



N 



(45) 



Furthermore, in the strong magnetic field limit, the Chern-Simons current projects to the 2d 
(x°, x^) plane as 



J, 



M 

cs 



MN A 

e An 



(46) 



Therefore, in the strong magnetic field limit, the theta term in (42) reduces to a mass term 



of the 2d gauge field, providing an explicit realization of the 2d Schwinger model [2T1 [7] , with 
the effective Maxwell equations in the reduced 2d {x^,x^) plane being written as 



J 



M 



□ + 



2^ 



A 



M 



(47) 



In physical terms, we see from (45) that a spatial inhomogeneity in the theta parameter 



corresponds to a non-zero A^{x^), which is a spatially inhomogeneous chemical potential. We 
expect this to produce charge separation. In terms of the Schwinger effect, this is illustrated 
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in Figure |2} depicting a spatial region of non-zero 9. At the edges, the gradient is non-zero, 
which produces a spatially inhomogeneous electric field, as shown in the Figure. This leads 
to an accumulation of opposite charges at the edges of the region of non-zero theta. On the 
other hand, a time-dependence in the theta parameter corresponds to a non-zero A^{x^), 
which acts as a chiral chemical potential in the LLL projection limit. For example, consider 
an electric field turning on smoothly at some early time, and turning off smoothly at some 
later time, as shown in Figure [s} Then the A^[x^) field is approximately constant over the 
time of nonzero electric field, and this is known to drive a current < > that is linear in 
A^[x^) [20], in agreement with the chiral magnetic effect relation (39). 



2 Chiral Magnetic Spiral 

In this Section we elaborate more on the dimensional reduction to two-dimensions due to 
strong magnetic field. In particular, we will show that the expectation values of vector and 
axial currents, < > and < >, can be expressed in terms of various fermionic bilinears 
whose dynamics are governed by a two-dimensional Lagrangian. The two dimensional axial 
anomaly reproduces the well known chiral magnetic effect. Furthermore, in addition to the 
chiral magnetic effect, the universal dynamics of two-dimensional chiral fermions implies the 
existence of additional currents which are transverse to the magnetic field, and have a spiral 
modulation along the direction of the field: the "chiral magnetic spiral" [22]. 



2.1 Basic setup and dimensional reduction 

We decompose the 4-component spinor in terms of eigenstates of the chiral projectors, Pr^l = 
^(1 ± 7^), the spin projectors P-^^^ = ^(1 ± S^), and the momentum direction projectors 



- i(l±7V) 

0^,3 _ fliag(^^3^ 



P. 

7S' 



The longitudinal spin operator is S" 



diag(cT , 0" ), and 



-0" ). We can write the 4-component spinor field as 



(4J 



The four-component spinor can be decomposed into two-component sub spinors in various 
ways. The chirality and spin decompositions respectively are 



R 



R. 



L 
L 

R. 



(49) 



where ± denotes the direction of motion along x^, the direction of the magnetic field. The 
corresponding four-dimensional currents can as well be decomposed in terms of chirality and 
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spin sub-spinors. The vector current \1/7^\1/ has the decomposition 

= ijjRipR- ipL ipL = -0tr^04 + 04r^0t 
j3 = ^nr^n + ijLr^L = h^'h + h^''h (so) 

The axial current = ^'-ff^'j^'^ has a similar form: 

J5 = i)RlpR + IpL IpL = + 

Here we have defined the two-dimensional gamma matrices as = = — za^, = a^. 

Let us consider a generic chiral Lagrangian 

where the interaction term does not couple left and right sectors. Since 7*^7^ is also block- 
diagonal in the chiral basis, the left and right sectors are completely decoupled and can 



Jl = V'Rr>i?-^Lr^^L = 0t0t-0>i (51) 



be treated as independent "fiavors" in the Lagrangian (52). Each chiral sector has its own 



associated current Jr l = 1/2(J^ ± J^) which are given by 

J] = aftpfijf 

4 = ^PfT^^Pf (53) 

Here f = L,R, and aR = —ai = 1. Even though the left and right currents seem to be 
completely independent, the conservation of vector current d^J^ = still holds. 

Now let us consider the lowest Landau level projection in the presence of a very strong 



magnetic field B directed along X3 direction as in Section lA If the magnitude of the 
magnetic field is the largest scale in the problem, then the transition from the lowest Landau 
orbit to an excited orbit will be exponentially suppressed and motion along the transverse 
X1X2 plane will be frozen. The kinetic term in the Lagrangian ([52|) then becomes: 



f=R,L f=R,L 

(54) 

by using the definition of two-dimensional gamma matrices above. This is the kinetic term 
of a two-dimensional Lagrangian. Let us further assume that the dimensionally reduced 
interaction term Cintjlipfizji)] is invariant under the two-dimensional chiral rotation: 

tfjfiz, t) e*^'^^("'*V/(^, t) (55) 
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for an arbitrary function (f{z,t). This two-dimensional notion of chirality generated by 
should not be confused with the four-dimensional chirality generated by 7^. The former acts 
on each right and left subspinor separately. Also, the four-dimensional chiral Lagrangian 



(52) does not have a term that couples right and left sectors so the system never develops 



a condensate such as < ipRipL >, and the four-dimensional chiral symmetry is never broken 
in our consideration. However, the dimensionally reduced system may exhibit dynamical 
breaking of the two-dimensional chiral symmetry (55)^ As an example let us consider the 



decomposition of the four-dimensional current-current interaction Cint 
(0, 3) components of the interaction become 



5^1- 



The 



(56) 



In two-dimensions this corresponds to a Thirring interaction. Similarly, the transverse com- 
ponents 



J}J 



n 



Jpf2 ^ i^f^Pf)' + ii^ftT.tff 



(57) 



generate a chiral Gross-Neveu or Nambu-Jona-Lasinio (NJL) interaction in two- dimensions 



Note that both (56) and (57) are invariant under (55) and exhibit chiral symmetry 



breaking in two-diemsions. The emergence of the Schwinger model from dimensional reduc- 



tion of four-dimensional Maxwell-Chern-Simons theory investigated in Section 1.5 is another 



example. To sum up, the four- dimensional currents (53) are governed by a two-dimensional 
chiral Lagrangian 



C2d= Yl V'/r^5M^/ + Ani,/[V^/(x*0] 
f=R,L 

after dimensional reduction due to the strong magnetic field. 



(58) 



We now show that the Lagrangian (58) has certain model-independent properties con- 
cerning the expectation values of the fermion bilinears < tpfTipf >. These are the building 



blocks of the expectation values of the four-dimensional currents (50, 51). Note from (50) 



that < Jj > and < Jj > are expressed in terms of two-dimensional densities and currents 
of ijjf, while the perpendicular components < Jj >,< Jj > are expressed in terms of two- 
dimensional scalar and pseudoscalar condensates of ipf. 



2.2 Life in two-dimensions 

Besides the magnetic field, the other necessary ingredient for chiral magnetic and/or sepa- 
ration effects is ambient charge density. In particular a nonzero baryon chemical potential 
/i leads to the chiral separation effect and a chiral (four-dimensional) chemical potential /i5 
leads to the chiral magnetic effect. To keep the discussion general, let us consider nonzero 
chemical potential for both right and left sectors separately 

fl'R = fJ' + 7^ , /iL = yU - yU5 7^ (59) 

^Spontaneous symmetry breaking in two dimensions is a delicate subject whose details are beyond the 
scope of our topic. It suffices to mention that in the limit of large number of flavors it can be realized, for 
example in the Gross-Neveu model [53]. 
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which can be reahzed by adding the term /i/V'j'?/'/ to the Lagrangian. In two dimen- 
sions the same term can be generated by a special local chiral transformation with a linear 
dependence on the spatial coordinate z in the exponent: 

^^ = e-^r'^/%. (60) 

To see this it is sufficient to observe 

^'f{tr'd,)^'f = ^f{tr'd,)^f + fif^}^f (61) 

where we used the gamma matrix identity r^F^ = r°. Since we assumed that the interaction 



term is invariant under any local chiral rotation of the form (55), the transformation (60) 
does not generate any other term than the coupling to the chemical potential. Note that this 
feature is special to two dimensions and does not generalize directly to higher dimensions. 
It is also crucial to have a chiral Lagrangian in four-dimensions. This would not hold for a 
massive fermion for instance. 

Let us start with the density < ip'jipf > and the reduced current along the z direction 
< ipfT^ipf > which constitute the four-dimensional density (< >) and current along the 



magnetic field (< Jj >). In Section (1.3) we have seen that the axial anomaly leads to an 
anomalous density in the presence of a chemical potential. In two-dimensional language, the 
axial anomaly 

Cm <Jf5> — omAjn 

< jf^ > = < ^^;fT''T'^Pf >= e^^^ < ^/L^^/^^ > (62) 

immediately reproduces < ^J-^/ >= yU//7r once eAj is identified with the chemical potential 
/i/. Alternatively the same result can be obtained directly from path integral by observing 
the transformation (60) applied to the renormalized charge density < i^'fi^'f >= creates the 
anomalous term fi j Jti 2A\ ■ 



ITT ^ ^ 

Once the projection to the lowest Landau level is implemented, only the spin up com- 
ponent 0-f- survives. Therefore for a right (left) handed spinor, only the positive (negative) 
momentum component contributes to the anomalous density, reducing its value by a half: 

< i^u'ipR > < R*,R+ > = 7;~77^ (LLL projection) 

zn zn 

<i){il)L> <L*_L_> = — ^ (LLL projection) (63) 

27r 27r 

The overall factor |^ is the density of the lowest Landau level in the transverse plane. There 
is no effect of the two-dimensional anomaly on the current ipfT^tpf in general. However, the 
lowest Landau level projection leads to a nonzero anomalous contribution to the current as 
well: 

<7/>«P^^> =<V't^a3V«> ^ <RlR+>=^^ 

ZH ZTl 

<^LT'i;L> =<ijlas^L> ^ -<L*_L_>=-^^ (64) 

Zn ZTT 
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The strong magnetic field selects a particular spin. Provided that there is some excess 
charge in the medium, the magnetic field induces a current for each chiral sector along its 
direction. When we transform the right and left currents back to vector and axial currents, 
we see the chiral magnetic and chiral separation effects in their conventional form: 

< J' >=< J° >= ^ < Jl >=< J° >= ^ (65) 

Now consider another set of bilinears, those forming the scalar < ipfipf > and pseudo- 
scalar < i/jfir^ijjf > condensates. They constitute the transverse components < > oi 
the four- dimensional currents. In the semiclassical limit where the number of fermion fiavors 
is large, the two-dimensional systems generically exhibit dynamical breaking of the chiral 



symmetry (55) and the system typically acquires a nonzero scalar condensate < ip'ip' >= 
m 7^ at zero chemical potential. Since this is also a mass term, it means there is a gap 
in the energy spectrum. The existence of the gap in the energy spectrum in one spatial 
dimension lowers the free energy at low temperatures as it "pushes" the Dirac sea further 
down in the energy spectrum. Once a finite chemical potential n is turned on, all the states 
with energy lower than fi will be occupied and the optimal configuration would be to open a 
gap right around to push the occupied states down and lower the free energy. This is the 
celebrated Peierls instability [25]. and has broad consequences in condensed matter physics. 



Our assumption of the invariance of the interaction under the chiral transformation (60) is 



sufficient to see that this scenario is indeed realized. The effect of (60) on the associated 
Dirac Hamiltonian is 

Htjj' = -iT'd,^/j' + Hint^' = {H- fi)^. (66) 

Therefore it is always possible to shift the energy spectrum, and hence the gap, by fi with a 
chiral rotation. This is the relativistic version of the Peierls instability and it is proven for 
instance explicitly in the 2d Nambu-Jona-Lasinio model [26], and is expected to be ubiquitous 
in two-dimensional systems with continuous chiral symmetry. 



The chiral transformation (60) mixes the scalar and pseudo-scalar condensates in the 
following way 

ip'fi^'f = cos{2 jj,fz)ijjfipf — sm(2 jj,z)ijjfiT^ipf (67) 

so the existence of a nonzero condensate < 'ip'f'ip'f >= rrif ^ at zero chemical potential 
generalizes into finite chemical potential as 

< ^p'^p >= mcos{2nz) , < ipiT^ip >= —msm(2ij,z) (68) 



since (67) must hold for any z. This modulated scalar/pesudo-scalar condensate is referred 
to as the "chiral spiral" [23 |2H] ■ In the 2d Nambu-Jona-Lasinio model the chiral spiral is 
indeed the thermodynamically preferred phase at low temperatures [2S]- The chiral spiral 



translates into transverse components of the four-dimensional currents modulated in z: 

< >= cr cos(2/i/jz + (pa) < JI >= -cl cos{2fiLZ + 4>l) 

< J\ >= Cr sm{2fiRZ + (Pr) < Jl >= Cr sin(2/ii2; + 0l) (69) 
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which we call the "chiral magnetic spiral" [22]. Here the amplitudes cr^l depend on the 
particular two-dimensional Lagrangian and are functions of the temperature, magnetic field 
and possibly other parameters in the model. However the chemical potential and space 
dependence of the currents is universal. (t)R^L are relative phases of the left and right chiral 
condensates. 

It should be emphasized that, as opposed to the longitudinal currents, the chiral magnetic 
spiral mixes up and down spin components. This can be seen in the spin decomposition of 
the currents (50, 51). In the lowest Landau projection, these pairings between spin up and 
down spinors correspond to excitations which can be described as a pairing of a particle with 
momentum /ij and a hole with momentum — yU/. This is because the particle and hole have 
opposite charges and therefore opposite spins in the lowest Landau level projection. Thus the 
excitation itself has momentum ±2/i/. This explains why the currents have the sinusoidal 
modulation in the z direction. 

In heavy ion collisions, the chiral magnetic spiral can induce both out-of-plane and in- 
plane fluctuating charge asymmetries (the explicit separation of out-of-plane and in-plane 
fluctuations has been performed [22] on the basis of STAR data [301 EI])- In the absence 
of topological fluctuations (/is = 0), at finite baryon density (/i ^ 0), and in the chirally 
broken phase, the charge current has only transverse components, and the charge asymmetry 
will fluctuate only in-plane. It should be kept in mind that the presence of magnetic field 
increases the chiral transition temperature [H]. If topological fluctuations are present in the 
chirally broken phase (e.g. due to the presence of meta-stable r]' domains [32J), the CME 
current can be carried by the chiral magnetic spiral. The chiral magnetic spiral has also been 
seen in a holographic study [33] in the framework of Sakai-Sugimoto model of holographic 
QCD 



3 Fermions in an instanton and magnetic field back- 
ground 

The previous discussion, along with many other papers, has presented the chiral magnetic 
effect as the flow of electrical charge as the result of some externally produced chirality 
imbalance, represented by a non-zero /is. In this Section we consider the situation in which 
this chirality imbalance is produced not by an explicit /is, or by a time-dependent A3, but 
instead as a result of a topologically non-trivial gauge background like an instanton. Since 
quarks carry both electric and color charge, they couple to both electromagnetic and gluonic 
gauge fields. In this Section we discuss some features of the spectral problem for fermions in 
the combined background field of a strong magnetic field and an instanton [35]. To illustrate 
the effect most clearly we take a single instanton in SU{2). We are motivated by situations 
in which quarks experience both types of fields, such as in dense astrophysical objects such 
as neutron stars and magnetars, and in heavy ion collisions such as those at RHIC and at 
CERN [SlEnilSZ]. 

We are also motivated by recent lattice QCD analyses [3H1 [39], SOl HH HSl |13] , which provide 
important numerical information about the Dirac spectrum in both QCD and magnetic field 
backgrounds. Analytically, while the effect of each individual background is very well known, 
their combined effect turns out to be quite intricate. In these lattice studies, certain matrix 
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elements associated with chiral effects receive dominant contributions from zero-modes and 
near-zero-modes, so we pay particular attention to the low end of the spectrum, and show 
that certain generic features have a very simple analytic explanation. 

As discussed already, a magnetic field introduces a Landau level structure to the fermion 
spectrum, in which the zero modes of the associated two-dimensional Euclidean Dirac oper- 
ator have definite spin, aligned along the magnetic field [HI |12]. For a constant magnetic 
field on a torus, as appropriate for lattice QCD analysis, this has been studied recently in 
im HSl HG] . The appropriate formalism is that of the magnetic translation group [H] . In a 
gluonic field with nontrivial topological charge (for example, an instanton), the fermion spec- 
trum of the four-dimensional Euclidean Dirac operator also has zero modes, with chiralities 
determined locally by the local topological charge of the gauge field |1H1 HHl EOl EH [52]. For 
a single instanton the fermion spectral problem has a conformal symmetry [53l and the 
zero modes are localized on the instanton, falling off as a power law with Euclidean distance. 
The conformal symmetry is broken by the introduction of a magnetic field, and now the zero 
modes develop an asymmetry, falling off in Gaussian form in the plane transverse to the B 
field, but as a power law in the other two directions. This basic asymmetry is an important 
feature of the phenomena of magnetic catalysis [9] and the chiral magnetic effect [H El El HI [5] , 
as sketched in Fig. |4} 




Figure 4: A sketch of the topological charge density, q oc trj^f^^^^^,, for a single instanton 
[red], and the density of the quark zero mode [grey]. On the left, there is a single instanton, 
and both densities fall off as power laws, with q falling off faster. On the right, with the 
introduction of a magnetic field, the topological charge density is unchanged but the zero 
mode density is distorted into an asymmetric shape, localized along the direction of the 
strong magnetic field. 



3.1 Euclidean Dirac Operator 

To discuss a combined background of an instanton and a static magnetic field we use Eu- 
clidean Dirac matrices, instead of the Minkowski ones used in the previous Sections. Our 
conventions follow those of [51], expressing the 4x4 Dirac matrices, 7^, for /i = 1, 2, 3,4, in 
terms of the 2x2 matrices = (1, and = (1, ia) = aj^, [here a are the usual 2x2 
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Pauli matrices]: 



Thus, the Euchdean Dirac operator can be expressed as 



where the covariant derivative, = — iA^, is written with a hermitean gauge field, 
Afj,, and is the Euchdean time coordinate. We write the gauge field as a sum of a 
non-abelian part, A^, and an abelian part, a^: 

A^, = A^ + (72) 

with the respective coupling constants absorbed into the gauge fields. The Dirac operator is 
anti-hermitean, so we write (with A real) 

z^V'A = AV^A (73) 



Since {75, _0} = 0, we can take A in (73) to be non-negative, with the negative eigenvalue 
solutions simply given by ip-x = '~f5'ip\. This means that we can effectively discuss the zero 
modes (A = 0) separately, and for the nonzero modes (A 7^ 0) we consider the squared 
operator: 

(z0)'^A= (^(f^ ^?^)^a = A2V'a (74) 

The positive chirality sector, x = +1; is described by the operator DD\ while the negative 
chirality sector, x = is described by the operator D'^D. We can write these operators as 

X = +l: DD^ = -Dl-^T,,a^, (75) 
X = -l: D^D = -Dl-^J^,,a^, (76) 

We have used [Dfj,,Di^] = — ij^ji,, where J-'^i, is the field strength associated with the gauge 
field A^, and the spin matrices a^i, and a^^ are defined as 



1 



a 



a^a„ - a^a^) , cr^z. = — (a^a^ - a^,a^) (77) 



In ( 75||76 ) we have used the properties [5Tj: a^aj, = + ia^y^ and a^a^ = + 



For non-zero modes [i.e., solutions to (74) with A 7^ 0], the operators DD"^ and D^D have 
identical spectra, for any background field. This is simply because we have an invertible 
map: suppose the 2-component spinor v satisfies D'^Dv = \^v. Then u = Dv is clearly 
an eigenfunction of the other operator, DD"^, with precisely the same eigenvalue: DD'^u = 
DD^Dv = \^u. Similarly, if u satisfies DD'^u = X'^u, then v = D'^u is an eigenstate of D'^D 
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with the same eigenvalue. Thus, when A 7^ 0, we can write the 4-component spinor solution 
in the form 



A 



^/'A = ( _i^t^^ ) where DD^ux = X^ux (7J 



or in the form 



^^A 

This is true for any background field: non-abelian, abelian, or both. 



= / A^^A j ^j^g^g D^Dvx = X\x (79) 



3.2 Magnetic field background 



For a constant (abelian) magnetic field, of strength B, pointing in the X3 direction, we have 
an abelian field strength fu = B, and so we find 



D^D = DD^ = -Dl - Ba3 



(80) 



where we have used the fact that o"i2 = ai2 = 0-3. Due to the subtraction term, —Bcx^, it is 
possible to have zero modes, and since DD'' = D'^ D these zero modes occur in each chiral 



sector. More explicitly, we can make a Bogomolnyi-style factorization similar to (12) and 
write 



Dl - Ba, 



-dl 



dl - {Dl T iD2) {Dl ± iD2) ± B 



-dl - dl - D^D, 



±B - Bas 



'1) 



For zero modes, we take c^s = ^4 = 0, and with i? > we choose the upper signs to ensure 
normalizable modes. For example, in the symmetric gauge where the abelian gauge field 



B 



-X2,Xi,0,0) 



(82) 



the zero modes can be expressed in terms of the normalizable solutions to {Di + iD2) u = 0: 



V^o = g{zi)e 






\V 



or 



ipo = g{zi)e 



-B\zi\y4 




1 



(83) 



Here g{zi) is a holomorphic function of the complex variable zi = {xi + 1x2). Both sets of 
zero modes have spin up, aligned along the B field; this is just the familiar lowest Landau 
level projection onto spin up states. Note also that the zero modes have the characteristic 
Gaussian factor in the {xi,X2) plane, transverse to the direction of the magnetic field. This 
factor is the origin of the distortion sketched in the right frame of Fig. |4j 

The number of zero modes per unit two-dimensional area [in the {xi,X2) plane] is given 
by the Landau degeneracy factor, the magnetic flux per unit area: B/{27i). In fact, even for 
an inhomogeneous magnetic field B{xi,X2), pointing in the 0:3 direction, the number of zero 
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modes [of each chirahty] is determined by the integer part of the magnetic flux (this is the 
essence of the Aharonov-Casher theorem [H]). For example, on a torus |12j : 

N+ = N^ = — I d^xB (84) 
27r J 

The higher Landau level states are the same for both spins, because {—D_D^ + B) and 
(— — B) have identical spectra, apart from the lowest level, which only has spin aligned 
along the magnetic field. 



3.3 Instanton background 

For an instanton field, A^, the (non-abelian) field strength F^y is self-dual [that is: F^y = F^^, 
where the dual tensor is defined: F^jy = ^efj,i,apFai3]. Then the anti-self-duality property of 
a^j, [that is: = -^^.upa^pa] implies: 

X = +l: DD^ = -Dl (85) 

X = -l: D^D = -Dl-^F,,a,, (86) 

Since —D^ is a positive operator, this means that for an instanton background there can be 
no zero mode in the positive chirality sector. On the other hand, due to the subtraction term, 
—F^iuCpu, in -D^-D, it is possible to have a zero eigenvalue solution in the negative chirality 
sector, and it has the form 



iIjo = , where Dv = (87) 

[For an anti-instanton, an anti-self-dual field with F^,^ = —F^^y, the zero mode lies in the 
positive chirality sector, because cr^j, is self-dual: a^i, = \^p,upa'7pa] For a general non-abelian 
gauge field A^, which is neither self-dual nor anti-self-dual, the Atiyah-Singer index theorem 
[561 152] states that the difference between the number of positive and negative chirality zero 
modes is given by the topological charge of the gauge field: 

N,-N_ = ^— f d^x F''.,.Ft 



327r^ 



Here we have written F^^ = F^^T"", with generators normalized as tr (T°T*) = For gauge 
group SU{2), with fermions in the defining representation, we take generators = ^r" in 
terms of the Pauli matrices r , and we can write the single instanton gauge field [55] , centered 
at the origin, in the regular gauge as 

K = (89) 

^ x^ + 

where p is the instanton scale parameter, and rj'^^ is the self-dual 't Hooft tensor [ISJ 1^ . 
The topological charge density is 



1 ~ 192 p^ 

^^^^ " 32^ ^'^"'^^'^"^ " (x2 + p2^4 (90) 
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There is a single zero mode [IHl HHl ISDl [iSl |52] , also localized at the origin, with density: 

64 p2 



(X2 + p2)3 



(91) 



These densities both fall off as power laws, with scale set by p, but the topological charge 
density is more localized, as indicated in the left-hand frame of Fig. |4| The nonzero modes 
are given by (78) or (79), and we note that the spectra are identical in each chiral sector, 
apart from the zero modes. 



3.4 Combined instanton and magnetic field background 

Physically, an instanton field projects the zero modes onto a definite chirality, while a constant 
magnetic field projects the zero modes onto definite spin, aligned along the direction of the 
magnetic field. When we combine the two background fields, both a non-abelian instanton 
field F^y and an abelian magnetic field /12 = 5, there is a competition between the two 
projection mechanisms, and the outcome depends on their relative magnitude, as we show 
below. Technically speaking, the instanton zero mode has a specific ansatz form that unifies 
space-time and color indices, while the magnetic zero modes have a natural holomorphic 
structure, and these two different ansatz forms do not match one another. The competition 
between these two ansatz forms makes the combined problem nontrivial. For an instanton 
field, since the field falls off as a power law, all eigenmodes also fall off with power law 
behavior. On the other hand, once a constant magnetic field is introduced, for example in 
the gauge ( [82| , all the eigenstates (even those in the higher Landau levels) have a Gaussian 
factor exp(— i?|2;ip/4) that localizes the modes near the axis of the magnetic field. This is 
the reason for the distorted density in the right-hand frame of Fig. 1. In the extreme strong 
magnetic field limit this leads to a dimensional reduction to motion along the magnetic field, 
with interesting physical consequences such as magnetic catalysis [9j and the chiral magnetic 
effect dllllE]. 

Concerning zero modes, we begin with a simple but important comment: in the index the- 



orem (88), the magnetic field makes no contribution, since with the field strength decomposed 
into its non-abelian and abelian parts, J^^^ = F^^, + f^^, we have 

+ (dim) f^J^, (92) 

= ii(F^,F^,) (93) 



where dim is the dimension of the Lie algebra representation of the non-abelian gauge fields. 
The cross terms vanish since the Lie algebra generators T" are traceless, and the f^uf\iu 
term vanishes since there is no abelian electric field. For example, if there is no nonabelian 
field, just an abelian magnetic field, then the topological charge clearly vanishes, and the 



index theorem (88) is consistent with the fact that DD'' = D'^ D for an abelian magnetic 



background (recall (80)), so that there are the same number of zero modes in each chiral 
sector. Now, with both background fields present, we find 

DD^ = -Dl-Bas (94) 



.2 1 



DW = -D'^~-F^,a^,~Bas (95) 
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Notice that the eigenvalues of DD'^ are simply those of the scalar operator — -D^, with a spin 
term ±5, as can be seen clearly in Figure 3^ The fact that there is a subtraction term from 
the positive operator — -D^ in both chirality sectors tells us that it is possible to have zero 
modes for each chirality, but their number will depend on the relative magnitude of F and B. 
In the next Section we study a specific model where we can quantify this precisely. Another 
important implication is that we may also have some "near-zero-modes", where the F and 
B subtractions do not exactly cancel the lowest eigenvalue of — -D^, but lower the eigenvalue 
of DD'^ or D'^D to near zero. 



3.5 Large instanton limit: Covariantly constant SU{2) instanton 
and constant abelian magnetic field 



In the very strong magnetic field limit, where the magnetic length, l/y/B, is small compared 
to the instanton size p, we expect a significant distortion of instanton modes and currents. 
In this limit we can make a simple approximation that reduces the problem to a completely 
soluble system. 

In the large instanton limit, we expand the instanton gauge field as: 

- J2^>- + • • ■ (96) 

To leading order in such a derivative expansion, the non-abelian gauge configuration A'^^{x) 
is self-dual and has covariantly constant field strength: F^^^ = —^Tj'^y In this limit we can 
make an SU{2) "color" rotation, along with a choice of Lorentz frame, to make the instanton 
field diagonal in the color space (we choose the direction), so that the field is self-dual, 
covariantly constant and quasi-abelian. Defining the instanton scale F = the combined 



gauge field, including also the abelian magnetic field as in (72), can be written as: 



F B 

A^, = - — (-X2, Xi, -X4, X3)r^ + y (-3^2, Xi, 0, 0)l2x2 (97) 

This gauge field is fully diagonal and moreover is linear in x^, so the problem is analytically 
soluble (this is the basic premise of the derivative expansion). The only nonzero entries of 
the field strength tensor are 



jr^2 = -Ft'^ + BI 



B-F 
B + F 



- -Fr^ - (98) 

In the absence of the magnetic field the field strength is self-dual, J^i2 = J^34, but a nonzero 
magnetic field breaks this symmetry. The topological charge density is (recall the normal- 
ization of the generators) 



1 ^„ ~, 4(2F)2 F 



'2 



32^2 i^^ 32;r2 27r2 
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A natural question to ask is: in such a constant field strength background, the wave 
functions have a Gaussian spatial dependence in the plane transverse to the direction of the 
field, characteristic of the Landau problem, so how can we recover the power-law dependence 
of the zero modes in an instanton background? This happens as follows. Recall that with 
the appropriate ansatz for the zero mode, the zero mode equation reduces to a first-order 
radial equation 



^0 



3r 



^0 



(100) 



l + r2'^" ' (l+r2)3/2 
where r is the Euclidean distance. At short and long distances, this zero mode behaves as 

3 



1 - -r^ + 
2 

+ ... 



r -> 



r — )■ oo 



(101) 
(102) 



On the other hand, if we make the above choice (97) of a constant field to represent a large 
instanton, we have instead the zero mode equation 



i/j'q = -3r Ipo 



^0 



-3r2/2 



(103) 



which has the correct short-distance behavior but which is Gaussian rather than power law 
at large distances. 




Figure 5: Succesive approximations in the derivative expansion of the instanton background 
interpolate between the exponential form of the zero mode [dashed line] to the power-law 
form [solid line] of the exact result. 

We can recover the correct power-law behavior by restoring the instanton size parameter, 
including the sub-leading terms in ( [96| , and expanding the zero mode equation as 

V'd = -3r (l - + - + - . . . ) t/^o (104) 



21 



The leading term represents the leading term of the derivative expansion of the instanton 
field, and is solved by the Gaussian factor. Absorbing this Gaussian factor by writing ipQ = 



-3rV2 



Xoi'f'), the resulting equation for Xo{f) is 



x'o = 3r (r^ - + - + . . . ) xo 



(105) 



which suggests writing xo 
derivative expansion, we obtain 

ipoir) = exp 



o(r). Continuing this process order by order in the 



So ^4 ^ fi 

ly I 7^ 7^ 

2 4 6 



exp 



ln(l + 



(106) 



thereby recovering the correct power-law decay from the derivative expansion of the instanton 



background, as illustrated in Figure 3.5 



3.6 Dirac spectrum in the strong magnetic field limit 

To study the Dirac spectrum with both a magnetic field and an instanton we consider the 
2x2 operators and DW in ^75f76^ . Notice first that 



J^^lu (Tf^u = 2 ( J'i2 - -^34) 0-3 (107) 

•^M-CT^. = 2 ( + 7-34) fX3 (108) 

It is convenient to factor the 4-dimensional Euclidean space and consider separately the 
{xi,X2) plane and the (^3,^4) plane. Then in the (xi,X2) plane we have a (relativistic) 
Landau level problem with effective field strength {B — F) in the = +1 sector, and with 
effective field strength {B + F) in the = —1 sector. In the {x3,X4) plane we also have 
a (relativistic) Landau level problem, now with effective field strength —F in the = +1 
sector, and with effective field strength F in the = —1 sector. In the {xi,X2) plane the 
sign of the effective field strength depends on which of 5 or F is larger, and in the strong B 
field limit, both B ± F are positive. 

When B > F, both [B — F) and {B + F) are positive. Thus, each color component of 
J^i2 is associated with a positive "magnetic" field. On the other hand, for J-34, the = +1 
sector has a negative field strength, while the = —1 sector has a positive field strength. 

We first consider the = +1 case. Then J^i2 = (B — F), T^i = —F, J^^,yO-^u = '^Ba^, and 
J^tyCT/iz/ = 2(5 — 2F)as. With a positive field strength the normalizable zero state is given 
by {Di + zDa) u = 0. But since ^^34 is negative, we factorize the corresponding covariant 
derivatives in the opposite order, in order to obtain a normalizable state annihilated by 
{D^ — iD^). Thus, we have, for chirality x = if? respectively: 

DD^ = -{Di- 1D2) {Di + tD2) - {Ds + iDi) {D3 - tD^) + B - Ba^ 

D^D = - (Di - iD2) (Z^i + tD2) - {D3 + iD^) (D3 - 2^4) + B-{B- 2F)a-i 

(109) 

This shows that there is a zero mode, when the spin term Ba^ cancels the B term from 
the Bogomolnyi factorization of the covariant derivative term. This occurs in the positive 
chirality sector, x = +1; and with spin up: (T3 = +1. 
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Now consider the = —1 case. Then J^i2 = {B + F), = F, Fp^uO'^y = 25(73, and 
FfjtuO'^v = 2(5 + 2F)(T3. All field strengths are positive, so we write, for chirality x = il; 
respectively: 



DD^ 
D^D 



{Di - iD2) (Di + iV2) - (^3 - iD^) (^3 + iD^) + {3 + 2F) - Ba^ 

{Di - iD2) {Di + iD2) - {Ds - iDi) (Dg + iDi) + (5 + 2F) - {B + 2F)<y^ 

(110) 



This shows that there is a zero mode, but now in the opposite chirality sector, x = and 
also with spin up: a^ = +1. 

To summarize: when B > F, the = +1 color sector has spin up zero modes with positive 
chirality, while the = —1 color sector has spin up zero modes with negative chirality. We 
can count the number of zero modes in each chirality sector by simply taking the product 
of the Landau degeneracy factors for the (xi,X2) and (x3,X4) planes, with the corresponding 
effective magnetic field strengths. Therefore, the corresponding Landau degeneracy factors 
give the zero-mode number densities (i.e, the number per unit volume): 



X = +1 : n+ = 
X = -l: n_ = 
The index (density) is the difference 



(B-F) F 

2tt 27r 
(B + F) F 

27r 2ti 



(ra = +1 , (T3 = +1) 
(rg = -1 , as = +1) 



(111) 
(112) 



EL 

'2^ 



(113) 



in agreement with the general index theorem (88), in view of (99). We also note that the 
total number density of zero modes 



n_|_ + n_ 



BF 
2^ 



(114) 



is linearly proportional to the magnetic field strength B. This is in agreement with numerical 
lattice gauge theory results [31]. 

It is worth emphasizing that if torus boundary conditions a;^ ~ + are imposed, 
then the fiuxes are quantized as BL^ = 27rM and FL^ = 211 N . Here M and N are positive 
integers. Consequently the index and total number of zero modes are given by also by integers 



index(^) = N+-N. = {N - M)M - {N + M)M = -2M^ 
total number of zero modes = (N + M)M + {N - M) = 2NM (115) 

Moreover, the instanton solution with constant field strength ( [96| is shown to be an exact 
solution of Yang-Mills equations on a four torus [STl EH]- Even with the inclusion of the 



magnetic field, (97) stays as a solution. Therefore on a on four-torus, there is no restriction 



on the magnitude of B and F (i.e. we do not have to consider a large instanton) and our 



counting of the zero modes (115) becomes exact for any integer M and N. 
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3.7 Physical picture: competition between spin and chirality pro- 
jection 

These results lead to the following simple physical picture. The instanton tries to generate 
a chirality imbalance but is neutral to the spin, whereas the magnetic field tries to generate 
a spin imbalance but does not affect the chirality. Depending on which is stronger, the zero 
modes have either a definite spin with a chirality imbalance {B > F), or a definite chirality 
with a spin imbalance [F > B). Also we see that in the former case, the total number of 
zero modes scales with B and is not equal to the index. 

More explicitly, for the B > F case, consider starting with just a strong magnetic field 
B, later turning on a weak instanton field. Without the instanton field, the zero modes and 



their degeneracy are given by the Aharonov-Casher theorem (84), so that the zero mode 
density is the Landau degeneracy factor B/{2tt) for each chirality sector. All the zero modes 
are spin up, as is familiar for the lowest Landau level (see Fig. |2]). There is an equal number 
of positive and negative chirality zero-modes, which is consistent with the index theorem, 
since the topological charge vanishes for a constant B field. Now consider turning on an 



instanton field F, with B > F > 0. We see from ( 111|112 ) that the effect of the instanton 
is to flip some of the chiralities: (^)^ positive chirality modes become negative chirality 
modes, leading to a chirality imbalance of in agreement with the index theorem (88). 
On the other hand, the total number of zero modes, grows linearly with the magnetic 
fleld when F is nonzero. 



3.8 Matrix Elements and Dipole Moments 

The computation of matrix elements is significantly simplified by using the Euclidean version 



of the Schur decomposition (18, 19) used earlier to relate the chiral magnetic effect to the 
two- and four-dimensional axial anomalies. Thus, introducing a small quark mass m, the 
propagator of the Dirac operator 'p + m is given by its Schur decomposition: 



1 Pit m 



(116) 



Note that DD'^ and D'^D have identical spectra, except for possible zero modes, so they can 
be viewed as square operators (matrices) of different dimension, as is clear when they are 
diagonalized in their respective eigenspaces. For a simple algebraic illustration, take D to be 
the 2-component column vector 

^=fil ^ DD'^={\ , Dt^ = 5 (117) 



2) \2 4^ 

Thus, DD'^ has eigenvalues and 5, while D'^ D clearly only has eigenvalue 5, the difference 
in rank being accounted for by the count of zero modes. 

The zero mode contribution to the propagator can be separated by writing it in one of 
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two ways, depending on which chirahty supports zero modes, 



-1 D 



1 _ / m^+DDt m2+DZ)t 



— ( (m m-^m2+DtD-^^) ^ m^+Dw\ (118) 

\ 1 nt m I ^ ^ 

An important set of quark bihnears involve the spin tensor E^j,: 



1 r n / a,„j 



S.. = ^l7..>l=(o- ,J (119) 

This representation makes clear the natural decomposition of E^uj, into its self-dual part {(Tf^i,) 
and its anti-self-dual part (o'^iy). The bilinears are 

< V^E^,^ >= tr (e^^^^^ (120) 

using the 2x2 sub-block structure of the propagator, it is straightforward to derive the 
relations: 

<^E,„^> = „tr,..(.--^+.--,^) (121) 
<#> = ™tr.,. (-.^ + -j^) (122) 
<^E„„^=^> = „t....(a-;;^j^-.-;^^j^) (123) 

<V.7=^> = „tr,,.(-j^--j^) (124) 

For apphcations to the chiral magnetic effect, we are interested in the magnetic and electric 
dipole moments: 

af = ^-eijk<i^^jki^> (125) 



af = <ilj^i4i^> (126) 

With a strong magnetic field in the direction, we concentrate on and a^, which require 
the spin tensors: 

Thus, 

^' m^ + DD^ ) + y^' m^ + DW ) ^^^^^ 

""' m^ + DD^ ) + '^^^^ \^' m-^DW ) ^^^^^ 
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The dominant contribution to the trace over the spectrum comes from the modes with low 
eigenvalues of DD'^ and D'^D. In the strong magnetic field limit, the zero modes and the 
near-zero-modes all have spin up, cjs = as expected. The dominant contribution to the 
electric and magnetic moments are therefore: 



2 \ / 2 

m \ I m 



m<^S,,^> - tr,,,( ^^^^^ l+tr...( ^^^^^ ) (130) 



2 \ / 2 

m \ I m 



m < ^E^^ij > ^ -tr2x2 , + ^^2x2 , n (131) 



m 



DD^ J """" \m^ + DW 



In particular, this means that in the strong magnetic field limit, we expect 

— ^ 1 , B ^ oo (132) 

<tpij > 

as has been confirmed in a lattice study 

For the magnetic dipole moment, the main contribution comes from the zero modes, so 
we simply count the degeneracies in the various sectors: 

Tv. , fB-F\ f F\ fB + F\ f F 

m < ^Si2^ > ^ — 7^ ■ 



27T J \2tI J \ 27T J \2'K^ 

which is linear in the magnetic field B. For the electric dipole moment, the near-zero-modes 
cancel, leaving just the zero mode contribution: 

™<..«.>.-(^^V^).(^)(£, 

(134) 



27r2 



which is independent of B, and negligible compared to BF, for B ^ F. [Note that (134) 
does not imply that there is a residual electric dipole moment when B vanishes, because (134) 
applies only in the B ^ F limit.] Thus, we see that the zero modes and near-zero-modes 
imply that 

< ^JjT.i2tp > (X. B , < tjjT.i2tJj > > < ^/'S34^/' > (135) 



This is in agreement with the lattice results of [o9] . We note that in a full QCD calculation 
with dynamical quarks there is an additional instanton measure factor that scales as m^^ , 
which should be take into account for these matrix elements. 

If we now consider the fluctuations in the electric dipole moment, we find a dependence 
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on B, because 



< ^'^34^ > = tr ( S34 Z.34 

IJ) + m ip + m 



1 „ 1 
rr 

_ f 1 .1 

- ^'^^^ \{m^ + DDrf + (m2 + DD^) ^^^^ ^^(m^ + DDt) ' + 
/ 1 . m2 1 

^ ^''^'{w^dW)^W^dw)) ^^^^^ 

where in the last step we have used the fact that the dominant contribution comes from zero 
modes and near-zero-modes, all of which have (T3 = +1. Thus, comparing with (133) we see 
that the fluctuation is linear in B 



< ^^34^ ^S34^ >^ ( ^^^^^ ) B (137) 



again in agreement with the lattice results of [39] 



4 Conclusions 

In this article we have presented several different perspectives of the chiral magnetic effect. 
The unifying theme is that the effect arises due to the spin projection of the lowest-Landau- 
level projection that occurs in a very strong magnetic field, in conjunction with a chirality 
projection that relates to the axial anomaly, or to the effect of a topologically non-trivial 
background field such as an instanton. The effect could occur in an abelian theory with an 
electric field parallel to the strong magnetic field. In this case, the roles of the charge and 
chiral chemical potentials are played by Aq and ^3, respectively, and the chiral magnetic effect 
is seen to be precisely equivalent to the 2d axial anomaly, which is itself the dimensionally 
reduced LLL projection of the 4d axial anomaly. Taking this dimensional reduction seriously, 
in a theory with a continuous chiral symmetry, we learn further from the 2d physics that 
there is a spiral condensate, an immediate consequence of the relativistic form of the Peierls 
instability. Finally, we considered the effect on light quarks of both a magnetic field and an 
instanton field, showing that the competition between the chiral projection in the instanton 
field and the spin projection in the magnetic field is responsible for the chiral magnetic 
effect. We demonstrated that this is consistent with the index theorem, and illustrated the 
mechanism with a soluble model by taking the leading derivative expansion form of the 
instanton field in which the magnetic length is much smaller than the instanton scale. It 
would be interesting to investigate the systematic corrections to this leading approximation 
using the Fock-Schwinger gauge |SSl EHl EI] representation of the instanton field. 

We thank T. Blum, D. Kharzeev and H-U. Yee for helpful discussions. This work was 
supported by the US Department of Energy under grants DE-FG02-92ER40716 (GD) and 
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